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Quasistatics is introduced so that it fits smoothly into the standard textbook presentation of elec- 
trodynamics. The usual path from statics to general electrodynamics is rather short and surprisingly 
simple. A closer look reveals however that it is not without confusing issues as has been illustrated 
by many contributions to this Journal. Quasistatic theory is conceptually useful by providing an 
intermediate level in between statics and the full set of Maxwell's equations. Quasistatics is easier 
than general electrodynamics and in some ways more similar to statics. It is however, in terms of 
interesting physics and important applications, far richer than statics. Quasistatics is much used 
in electromagnetic modeling, an activity that today is possible on a PC and which also has great 
pedagogical potential. The use of electromagnetic simulations in teaching gives additional support 
for the importance of quasistatics. This activity may also motivate some change of focus in the 
presentation of basic electrodynamics. 



I. INTRODUCTION 

Applications of electrodynamics may be in the static, 
quasistatic or general high frequency regime. Quasistat- 
ics is more or less neglected in most textbooks and the 
purpose of this paper is to present material that can fill 
this gap in a course on the level of for example Griffith's 
textbook.^ There are good reasons to do so including: 

1. The step from statics to general electrodynamics is 
huge in terms of its physical content. Quasistatic 
models are useful by providing intermediate levels 
in the theory. Thereby several confusing issues that 
result from the very condensed standard deriva- 
tion of Maxwell's equations may be avoided. Ques- 
tions concerning Coulomb's and Biot-Savart's laws 
in non-static situations should be addressedi2i2i4 
So should also the appearance of dynamic electric 
fields in regions where the magnetic field seems to 
be absent (like outside a long solenoid or toroidal 
coil with a time- varying current) i^iSiSiS Another 
confusing issue follows from the standard textbook 
motivation for the displacement current where it is 
introduced in order to make the equations of elec- 
trodynamics consistent with charge conservation. 
The physical impact of this term is remarkable and 
includes in particular electromagnetic waves in free 
space. But why should waves in free space without 
charges and currents follow from charge conserva- 
tion?— Quasistatics is useful in the discussion of 
these and many other questions. 

2. There are plenty of interesting phenomena within 
the quasistatic regime. Using the full set of 
Maxwell's equations are many times unnecessarily 
complicated since these equations can describe the 
most intricate electromagnetic wave phenomena in- 
volving short time-scales or high frequency. Such 
an analysis may be difficult and is not necessary in 
quasistatic situations. 

3. Many real world applications involve the numeri- 



cal solution of the time-dependent Maxwell equa- 
tions in three dimensions. However, quite often one 
is interested in phenomena where some quasistatic 
analysis is sufficient. This amounts to the replace- 
ment of an hyperbolic model with an elliptic (or 
parabolic) one which can be solved in more eco- 
nomic waysii2iiii4 

4. Basic quasistatic theory may be introduced in an el- 
ementary and simple way suitable for basic courses 
in electrodynamics. Like static theory also qua- 
sistatics may be given two equivalent formula- 
tions. The first is in terms of force laws, includ- 
ing straightforward generalizations of the static 
Coulomb and Biot-Savart laws. The second for- 
mulation is in terms of simple approximations in 
Maxwell's differential equations. 

This paper is organized as follows. Basic theory of 
quasistatics is presented in Section II. Subsections A and 
B contain the two equivalent formulations of quasistat- 
ics, one in terms of force laws and the other in terms 
of approximations in Maxwell's differential equations. A 
salient feature of quasistatics is instantaneous interac- 
tion at a distance. In subsection C the corresponding 
c — > oo limit of general electrodynamics is used to derive 
quasistatics. Subsections D, E and F include quasistatic 
theory for the electromagnetic potentials, the fields from 
a moving point charge and the quasistatic Poynting the- 
orems. Alternatives to the standard textbook derivation 
of Maxwell's equations are discussed in Section III from 
a quasistatic perspective. 

Sections IV and V concern some confusing issues that 
have been discussed in many contributions to this Jour- 
nal. In Section IV we consider the laws of Coulomb and 
Biot-Savart within general time-dependent theory.^ Suf- 
ficient and necessary conditions for the exact validity of 
these laws are given. Section V contains comments on 
the old question #6 of this Journal which concerns the 
fields outside a solenoid with a time-varying current^ 

Quasistatics is useful not only for providing improved 
understanding of basic theory but it is also important 
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for applications, in particular for numerical simulations. 
In Section VI we discuss the possibility of using elec- 
tromagnetic simulations in basic courses and exemplify 
with quasistatic equations for eddy currents. A couple of 
PDE-solvers which have been used for electromagnetic 
modeling on a PC is mentioned. A summary is given in 
Section VII. 



II. QUASISTATICS 

Quasistatics within electrodynamics refers to a regime 
where "the system is small compared with the electro- 
magnetic wavelength associated with the dominant time- 
scale of the problem" The fields are then propagated 
instantaneously so we are dealing with a kind of c — *■ (X) 
limit. We consider in this paper three major quasistatic 
models. These are EQS (electroquasistatics)ii, MQS 
(magnetoquasistatics)^^ and the Darwin modeliSiiiii. 
EQS includes capacitive but not inductive effects, MQS 
includes inductive but not capacitive effects while the 
Darwin model includes both capacitive and inductive ef- 
fects. The Biot-Savart law is valid in all three models 
while the Coulomb law is valid only in EQS. In MQS 
and Darwin there is also, besides the Coulomb field, an 
additional contribution to the electric field due to mag- 
netic induction. The source of this electric field is thus 
(9B /dt in Faraday's lawii& 



A. From the laws of Coulomb, Biot-Savart and 
Faraday to the EQS, MQS and Darwin models. 

Dynamical systems that proceed from one state to an- 
other as though they are static (at each fixed time) are 
commonly said to be quasistatic. For electromagnetism 
the static theory builds on the Coulomb and Biot-Savart 
laws together with the static continuity equation. Qua- 
sistatics would then simply be obtained by the allowance 
for time-dependence in the first two otherwise unchanged 
laws 



E(r,t) = 



p(r',i)R 



dr' 



B(r,t) 



Mo 



J(r',i) X R 



dr' 



where we use the notationaii 

R r - r', i? = |R| , R = R/i? 



(1) 



(2) 



(3) 



It would however appear strange to keep the static conti- 
nuity equation unchanged in this time-dependent situa- 
tion so we replace it with the usual continuity equation^^ 



dp 

dt 



V • J = 0, 



(4) 



The three equations (P), 10) and l@J constitute the EQS 
model that includes capacitive but not inductive effects. 
Charge may be accumulated in this model but this re- 
quire work and energy may as usual be associated with 
the electric field. Magnetic energy is however outside the 
scope of EQS because there is no magnetically induced 
electric field and accordingly no back emf. No work is 
then required to create a magnetic field by starting an 
electric current. 

We now like to include electromagnetic induction. 
Then also d'B/dt acts as a source of electric fields and 
the total electric field becomes the sum of two parts 



The first term is the Coulomb electric field 



Ec (r,t) 



1 



47reo 



(5) 



(6) 



and the second term is the Faraday electric field which 
may be defined by a Biot-Savart like integral expression 
as (cf. equation 



EHr,t) = --l 



OB {r',t) R , 



We note that Ep solves the equations 

dB 

V X Ef = — 7^, V • Ef = 
ot 

while the corresponding equations for Ec is 

V X Ec = 0, V • Ec = — 

eo 



(7) 



(8) 



(9) 



The expressions and {Tj) are the unique solutions of 
^ and 10 provided appropriate boundary conditions at 
infinity are usedi^ 

The Darwin model may be defined by the equations 
and (©-(EI). This is a quasistatic model that includes 
both capacitive and inductive phenomena. Note that for 
given current and charge densities we directly obtain the 
electromagnetic fields in terms of integrals without the 
appearance of any time-retardation. Instead the integral 
expressions in the laws of Coulomb and Biot-Savart play 
an important role in this dynamical model. 

The MQS model is obtained from Darwin if the usual 
continuity equation is replaced by the static one 



V • J = 



(10) 



Thus the MQS model may be defined by the equations 
@, ©, ®, O and It is different from both EQS 

and Darwin by including inductive but not capacitive 
effects. A confusing feature of MQS is that the very fun- 
damental continuity equation may be violated by MQS- 
solutions. Only stationary currents are allowed in MQS 
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and these cannot explain changes in the charge density, 
thus one cannot interpret the currents in MQS in terms of 
charge transportation. Of course, only sufficiently good 
approximations of solutions to the Maxwell equations are 
interesting in the real world so the continuity equation 
must still be almost true in some sense. Ampere's law, 
which imply equation (|10|l . is valid in MQS but not in 
EQS or Darwin (see the next subsection). Thus Ampere's 
law is not always valid in quasistatics. Griffiths and 
Heald remark that "The application of Ampere's law 
in quasistatic situations can be an extremely delicate 
matter" ? 

The textbook of Haus and Melcher— builds up an un- 
derstanding of electrodynamics by using both EQS and 
MQS. This is of particular significance in the relation 
between electromagnetic field theory and circuit theory. 
Then EQS involves capacitance features and MQS induc- 
tance features. For systems involving capacitance and 
inductance both models are needed. For such applica- 
tions it is crucial that capacitive and inductive aspects 
are not both important in the same spatial placcpSi The 
use of two complementary quasistatic models in the same 
physical system is clearly a complicating feature if we like 
to model the whole system numerically. We would then 
have to divide the whole spatial region into EQS and 
MQS subregions with appropriate continuity conditions 
at the interfaces. A better alternative may be to use the 
Darwin model which embrace all the physics contained 
in EQS and MQS, still being quasistatic. 

B. From Maxwell's equations to EQS, MQS and 
Darwin. 

Let us now consider the formulation of quasistatics in 
terms of differential equations. The starting point is gen- 
eral electrodynamics with the Maxwell equations: 

V-E=(l/eo)p (11) 



V-B = (12) 



TABLE L This table give two equivalent definitions of each 
of the quasistatic models EQS, MQS and Darwin 



Model 


Def.A, eq.# 


Def.B, eq.# 


EQS 




ijTT|i.iP|i.iPji.ijT^ 


MQS 


iPi.pi.p.i|7ji.iP|i 


ijTT|i.ip|i.in^.inii|i 


Darwin 




i(nj, 111211,10,1113, © 



and B-fields will be used in the present paper but it is 
straightforward to modify the expressions so that a for- 
mulation with D- and H-fields is obtained. 

Let us now formulate the quasistatic models EQS, 
MQS and Darwin as approximations of Maxwell's equa- 
tions Hll|l - H14() . The relations to the integral formulas 
in subsection A above will also be considered. The EQS 
model is obtained from the Maxwell equations simply by 
neglecting 9B /dt in Faraday's law while MQS is obtained 
by instead neglecting dE/dt in the Ampere-Maxwell law. 
Thus in EQS we have instead of (|13|) the equation 

VxE = (15) 

and in MQS the equation fl^ is replaced by the usual 
Ampere law2^ 

V X B = ^oJ (16) 

To get the Darwin model we do not neglect all of dE/dt 
but keep the Coulomb part of the E-field (defined by 
equation © above) and replaces the Ampere-Maxwell 
law with the Ampere-Darwin equationS^ 

VxB = ^o(^J + eo^) (17) 

Note an important difference between these two last 
equations. The Ampere law (|16|l implies equation I|1U|) 
and may violate charge conservation in dynamic situa- 
tions. The Ampere-Darwin equation, like the Ampere- 
Maxwell equation, is consistent with the continuity equa- 
tion Q. 



V X E 



V X B = ^0 J 



dB 

"dt 



(13) 



(14) 



These equations are complete and general as they stand. 
However, in the presence of polarizable/magnetizable 
media it is in practice very convenient to write them in 
a different way by introducing the D- and the H-fields. 
Then only the free charge and current densities appear 
explicitly in the equations. For formal simplicity just E- 



The electrostatic models may thus be defined with fo- 
cus on the laws of Coulomb and Biot-Savart (like in the 
previous subsection) or alternatively as approximations 
of Maxwell's equations. These definitions are summa- 
rized in Table The proof that the A and B columns 
define the same models involves only standard proce- 
dures found in basic textbooks. In electrostatics one start 
from the Coulomb law and obtain the divergence and 
curl of the electric field. By the Helmholtz theorem we 
have the equivalence (assuming appropriate conditions at 
infinity^") 



^ ^ and (nni 



(18) 
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TABLE II: The table includes several familiar laws from elec- 
tromagnetics as well as the Maxwell-Darwin equation and 
state if they are valid within the EQS, MQS and Darwin mod- 
els. 



Eq.# 


Equation 


EQS 


MQS 


Darwin 


CD 


Coulomb's law 


yes 


no 


no 




Biot-Savart's law 


yes 


yes 


yes 




Continuity eq. 


yes 


no 


yes 




V • J = 


no 


yes 


no 




Gauss law 


yes 


yes 


yes 




V B = 


yes 


yes 


yes 




Faraday's law 


no 


yes 


yes 


m 


Ampere-Maxwell 


yes 


no 


no 




V X E = 


yes 


no 


no 


m 


Ampere's law 


no 


yes 


no 




Ampere-Darwin 


yes 


no 


yes 



In the same way we of course have (^l © • In magneto- 
statics one starts from the law of Biot-Savart and obtain 
the divergence and curl of the magnetic field. One finds 
the equivalence 



© and (Cni) ^ (O and (0 



(19) 



The above relations trivially remain valid if we allow for 
time-dependence where time appears only as a parame- 
ter. However, in a time-dependent situation it is logical 
to use the general continuity equation Q instead of the 
static one (|10|) . Then instead of (|19|l we find 



lEl and gl) lO and (O 



(20) 



Thus the Ampere-Darwin equation is a mathematical im- 
plication of Biot-Savart's law combined with the continu- 
ity equation. The relation (|20|l is not included in most 
standard textbooks but has (more or less explicitly) ap- 
peared in several contributions to this Journali2iSSiS&. 

The equation Q for the Faraday electric field is for- 
mally analogous to the Biot-Savart law l^l for the mag- 
netic field. Formally similar to H19|) is the equivalence 



O and 



© and lO 



(21) 



Table |TT| includes some laws of electromagnetics and 
states if they are valid within EQS, MQS or Darwin. All 
the included equations are familiar from standard electro- 
magnetics with the exception of Ampere-Darwin's equa- 
tion ifT7|) . 



C. A limit of instantaneous propagation in 
Maxwell's equations 

In subsection A above we obtained the Darwin model 
from EQS by including Faraday's law in a simple and 



straightforward way. An alternative procedure is to 
follow the standard textbook derivation of Maxwell's 
equationsi but in the final step, when the displacement 
current is added to Ampere's law for consistency with 
the continuity equation, find the Ampere-Darwin equa- 
tion H17(l rather than the Ampere-Maxwell equation H14|l . 
This replaces the usual maximal assumption in the cor- 
recting termlS with a kind of minimal assumption. Both 
these derivations of Darwin's model may however seem 
rather superficial and arbitrary. It would be nice to obtain 
it from Maxwell's equations by using some more system- 
atic method. A basic feature of quasistatic approxima- 
tions to Maxwell's equations is the instantaneous prop- 
agation of fields. Thus it should be possible to consider 
quasistatics as some limit c ^ oo of Maxwell's equations. 
This is however quite a singular limit and the procedure 
must be further specified. We take the absence of time- 
retardation as being the most salient property of qua- 
sistatics. Let us express Maxwell's equations in terms of 
the potentials {V, A) so that the retarded time appears 
explicitly. We choose to use the Coulomb gauge 



V • A = 



(22) 



and follow the derivation in Nielson and Lewisi^ (the rea- 
son for not using the Lorentz gauge is considered soon). 
The equations for the potentials become (see Griffithsi 
p. 421), 



So 



w2 . S^A 



-fioJ + MoeoV — 
ot 



(23) 



(24) 



Time-retardation appears explicitly if we solve equation 
(|24|l for the vector potential in terms of an integral in the 
usual way. The omission of retarded time in this integral 
is the same as excluding the second order time-derivative 
so that 



V^A = -^oJ 



„dV 
Aio^oV — 
at 



(25) 



The model so obtained consists of the equations H22|l . 
(I23|l and (|25|l . These constitute the Darwin model in 
terms of potentials. The usual Darwin equations Hll|l . 
(221 : 03 and ^\ is obtained by use of 



E = -VV - 



dA 

~dt' 



B = V X A, Ec = -VF (26) 



Why didn't we use the Lorentz gauge? It is straight- 
forward to make the corresponding calculations also in 
that case. We then get a model approximating Maxwell's 
equations which seems pretty close to the Darwin model 
with the Ampere-Darwin law as one of the equations. 
However, Gauss law (|llf) is not obtained but a new term 
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appears in the corresponding equation. Certainly we pre- 
fer not to change Gauss law and this problem, as we have 
seen, does not appear if we use the Coulomb gauge in the 
approximation procedure. Using the Coulomb gauge in 
the Darwin model also have the nice consequence that 
the Coulomb part of the E-field is given by the scalar po- 
tential Ec = — Vy and the Faraday part by the vector 
potential Ep = —dA/dt. 



D. The potential representations of EQS, MQS and 
Darwin 

We use the Coulomb gauge. The magnetic field is writ- 
ten B = V X A for all models. For MQS the vector 
potential then satisfy 



(27) 



1. The force between two point charges depends not 
only on their separation but also on both their ve- 
locities and accelerations 

2. Furthermore, it is the position, velocity and accel- 
eration at retarded time of the other particle that 
matters. 

The second point is far more challenging than the 
first one. However, time-retardation vanishes within qua- 
sistatics and the approach with interacting point charges 
become quite simple and instructive. In the point charge 
approach to electrodynamic models wc automatically in- 
clude the continuity equation Therefore this section 
concerns EQS and Darwin but not MQS. Let us start 
with the EQS model. The fields from a point charge Q 
within EQS is 



or in integral form 



E(r,i) 



R 



(33) 



4tt R 



(28) 



For EQS and Darwin the corresponding equations are 
(113 or 

It takes some manipulation involving partial integration 
and the continuity equation to derive H29|l from (|25|) . The 
scalar potential is just the Coulomb potential for EQS, 
MQS and Darwin, 



V 



47reo J J J R 



-dr' 



(30) 



The electric field for EQS is conservative and is just the 
Coulomb field 



B(r,i) 



^oQ V X R 



(34) 



where the notations in J^Jl are used. The position of Q 
is r' which is a function of time and v ~ dr' /dt. From 
these equations we obtain Q and ||2J| while the continuity 
equation is implied by the point charge description. 

Let us now consider the Darwin model. Then equation 
(|34|l remains valid while l|33|l only gives the Coulomb part 
of the electric field 



Ec (r,t) 



Q R 

47reo i?2 



(35) 



The magnetically induced part 'Ep must now also be in- 
cluded. However, it is then convenient to start all over 
again using the electromagnetic potentials. The poten- 
tials are 



(31) 



while the electric field for MQS and Darwin also include 
the magnetically induced electric field so that 



E = -VV 



dA 
dt 



(32) 



E. Fields of a moving point charge 

We will now consider the quasistatic fields from a point 
charge in general motion. Coulomb's law is often taken 
as a starting point for electrostatics. It is conceptually 
simple to build on interactions between point charges in 
this way, so why not develop all of electrodynamics by 
generalizing this approach? This is discussed by Griffiths 
(chapter 2).-'^ Two major problems are 



V{r,t) 



AttsoR 



A{r,t) 



i^oQ V + V • RR 



47r 



2R 



(36) 



(37) 



The scalar potential H3()|l is just the Coulomb potential 
and the vector potential H37() may be found from 129|l . 
An alternative derivation is by the ansatz 



(38) 



where the scalar field 4> is determined by the Coulomb 
gauge condition l|22(l . It takes a little algebra to find 



' Stt 



V R 
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and thus (|37(l The electromagnetic fields is related to 
the potentials in the usual way (|26|l . The magnetically 
induced part of the electric field is 



E, 



-dA/dt 



(39) 



This expression will obviously involve the acceleration 
a ~ dv/dt oi Q. The Darwin model should thus follow 
from the potentials and (|37|l of a point charge and 
it is sufficient to check ©-(Ell- Only equation Q is not 
obvious and the easiest approach is to use the equiva- 
lent differential equations But these equations are 
trivially satisfied by using the Coulomb gauge and 
B = V X A. 

The standard textbook approach to electrostatics be- 
gins with Coulomb's law (|33|1 for a charge at rest. Some- 
times the corresponding approach to magnetostatics is 
used by taking as a starting point but now, of course, 
with a moving charge. This may in principle seem wrong 
since this is a time-dependent system. Most textbooks 
take this seriously and starts instead from the Biot-Savart 
law with stationary current. However, leaving the point 
particle point of view also makes magnetostatics more 
difficult than electrostatics. It comes as a relief that the 
magnetostatic differential equations (|12|l and (|16|l are for- 
mally similar to, and not much more difficult than, the 
electrostatic equations and (|15|) . Quasistatics is eas- 
ier than statics in this respect; it allows for the use of 

A conventional Lagrangian description of the electro- 
magnetic interaction between two or more charged par- 
ticles is possible only if time-retardation may be ne- 
glected. It is straightforward to give such a formulation 
constructed from the potentials and (|37|1 .'^" This re- 
sults in the so called Darwin Lagrangian that was first 
obtained by Oliver Heaviside in 1891i^ii2^ It was found 
again by C.G. Darwin in 1920 by an expansion of the 
Lienard-Wiechert potentials.^^ Rather surprisingly it 
then turns out that that terms up to the order (w/c)^ 
inclusive are to be kept in the Darwin Lagrangian (see 
Jacksoni^p.596). 



F. The Poynting theorem 

Let us follow the standard derivation of Poynting's the- 
orem while using the quasistatic models. I.e. we add the 
two equations obtained by taking the scalar product of 
the curl E equation with the B-field and the scalar prod- 
uct of the curl B equation with the E-field. We the get 
(after the usual manipulations) for EQS 



|(ieoi.^) + V.(i-ExB) 



for MQS we obtain 



ot 2/io A^o 



E J 



E J 



(40) 



(41) 



and for the Darwin model 

|,l..£g, + _l_B=,+v.(J-ExB + .„ff).-E^J 

(42) 

In EQS only the electric field is associated with energy. 
Building up a magnetic field costs no energy in this model 
due to the absence of a counteracting induced E-field. In 
MQS only the magnetic field has energy. Changes of the 
electric field is associated with changes in the charge den- 
sity. However, the continuity equation is not satisfied and 
there is no charge transport for which we could calculate 
the required energy. Finally the Darwin model includes 
both magnetic and electric energy. Note that only the 
Coulomb-part of the E-field is associated with energy. 
An important qualitative difference between EQS and 
MQS on the one side and Darwin's model on the other is 
the possibility of natural resonances in the latter. Since 
Darwin includes both capacitive and inductive features 
we may in principle use these equations to model, for ex- 
ample, some field theoretic manifestation of a LC circuit. 



III. IF MAXWELL HAD WORKED IN 
BETWEEN AMPERE AND FARADAY 

We consider in this section various procedures to find 
the full set of Maxwell equations starting from the laws of 
Coulomb and Biot-Savart. The quasistatic perspective is 
useful and the models EQS, MQS and Darwin will appear 
as intermediate stages. The title of this section refers to 
the possibility of introducing the displacement current 
before Farday's lawi^ 

Let us start with the standard textbook procedure of 
finding Maxwell's equationsi From the static laws of 
Coulomb and Biot-Savart we find, by use of the static 
continuity equation, the static limit of Maxwell's equa- 
tions. By introducing magnetic induction we then obtain 
MQS (i.e. " Electrodynamics before Maxwell" in Griffiths 
textbook). The final step, motivated by the need for con- 
sistency with charge conservation Q , is to introduce the 
displacement current. At this point a question mentioned 
already in the introduction may surface: Maxwell's equa- 
tions includes waves in free space but why should these 
follow from charge conservation! Quasistatics is useful 
for discussing this issue. It is indeed not necessary to 
introduce all of the displacement current to save charge 
conservation. It is sufficient to include the Coulomb part 
of itiS and then we get the Darwin model. At this point 
no surprising new physics appears. The step from Dar- 
win to the Maxwell's equations may then, to begin with, 
be motivated by symmetry and beauty of equations. The 
explicit appearance of the Coulomb electric field in the 
Darwin model is a rather unsatisfactory feature and it 
is so easily fixed by just replacing it with the total elec- 
tric field. Amazing new physics now appears and this is 
accompanied with beautiful new mathematical structure 
like the Lorentz invariance of both the Maxwell equations 
and of the trajectories of test charges. This is of course 
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surprising but, at least, the previous rather mystifying 
motivation for it has now been removed. 

Jammer and Stachel discuss what might have hap- 
pened if MaxweU had worked in between Ampere and 
Faraday.3^ As in the standard textbook derivation of 
MaxweU's equations one may first derive the static hmit 
of Maxwell's equations from the laws of Coulomb and 
Biot-Savart combined with the static continuity equa- 
tion. At this stage Maxwell might have added the dis- 
placement current before the discovery of Faraday's law. 
This would result in the EQS model which, according to 
Jammer and Stachel, is exactly Galilei invariant. They 
suggest in the article abstract that the discovery of Fara- 
day "would have confronted physicists with the dilemma: 
give up the Galilean relativity principle for electromag- 
netism (ether hypothesis), or modify it (special relativ- 
ity). This suggest a new pedagogical approach to elec- 
tromagnetic theory, in which the displacement current 
and the Galilean relativity principle are introduced be- 
fore the induction term is discussed." This approach is 
however less striking than it first seems. We like to in- 
clude the appropriate invariance structure not only for 
the field equations but also for the trajectories of test 
charges. Here one get problems with the Galilei invari- 
ance of EQS while the Lorentz invariance of Maxwell's 
equations is perfect. The problem of finding satisfac- 
tory quasistatic and Galilei invariant approximations to 
Maxwell's equations is not an easy onei^ 

A third approach to Maxwell equations is of more qua- 
sistatic nature. We avoid in this derivation the static 
limit of Maxwell's equation by allowing for (trivial) time- 
dependence in the laws of Coulomb and Biot-Savart 
(121). We obtain the differential equations of EQS (i.e. the 
equations 1)11(1 . H12|l . (|14|l ) and l|15|l by pure mathematics 
if we also use the (time-dependent) continuity equation 
Q. Maxwell's equations now follow directly when the 
Faraday law is introduced. Also in this case (like in the 
standard textbook derivation) a question may appear. 
Why do the magnetic induction experiments of Faraday 
imply electromagnetic waves and time-retardation? The 
answer is somewhat more hidden in this case. The elec- 
tric field in EQS is a pure Coulomb field Ec but when 
we include magnetic induction we also get the Faraday 
electric field Ep. It is now not clear if we should have 
E or Ec in (|14|l . In the latter case (the minimal as- 
sumption) we get the Darwin model and no qualitatively 
new physics appears (only what is needed to explain the 
experiments of Faraday). The step from Ampere-Darwin 
to the Ampere-Maxwell law has already been discussed. 

IV. THE LAWS OF COULOMB AND 
BIOT-SAVART IN TIME-DEPENDENT THEORY 

The laws of Coulomb and Biot-Savart provide a start- 
ing point for the static theory of electromagnetics. Let 
us consider the following two questions: 

1. How should the static laws of Coulomb and Biot 



-Savart be generalized to time-dependent theory? 

2. Consider the trivial generalization given by the 
equations ^ and (0) where time is included as a 
parameter. What is the significance of these for- 
mulas? 

These questions are usually not addressed in the text- 
books. The first question is considered by Jefimenko and 
the resulting formulas involves (as must be expected) in- 
tegrals where the retarded time appearsi^ Jefimenko's 
generalized laws of Coulomb and Biot-Savart are used by 
Griffiths and Heald to address the second question. They 
find thali^ 

(a) The generalized Biot-Savart law reduces to the 
standard one if d^J / dt^ = (then, from the conti- 
nuity equation, also d^p/dt^ — 0). 

(b) The generalized Coulomb law reduces to the stan- 
dard one if dJ/dt — (and thus also d^p/dt^ = 0). 

(c) The law of Ampere holds if d'^3/dt^ = and 
dp/dt = 0. 

These results are all sufficient conditions for the laws 
but only (c) is also necessary. Below necessary and suf- 
ficient conditions are formulated. It is natural to con- 
sider the validity of Coulomb and Biot-Savart in time- 
dependent theory from the perspective of quasistatics. 
We may divide electromagnetics into three regimes: stat- 
ics, quasistatics and high frequency. However, statics is 
just a particular case of the general Maxwell equations 
while quasistatics is an approximation. This motivate a 
third question : 

3. When does it happen that a solution in quasistatics 
is an exact solution to Maxwell's equations? 

It will be shown that a solution of the Darwin model 
also solves Maxwell's equations if and only if the current 
has the form 

J (r, t) = a (r)< + b (r) - e, (43) 

at 

where the Coulomb field is defined by JHl and the vec- 
torfields a and b satisfy V • a = V • b = 0. Formally we 
write this 

MaxweU -I- Biot-Savart Darwin + Eq.lO (44) 

Let us compare this with (a) above. The last term in 
(|43|l may appear new and unfamiliar in the context of 
being a part of the true current (this term but with op- 
posite sign is familiar as a part of the displacement cur- 
rent). However, this term is just the solenoidal part of 
the current. Thus consider any current density J and 
write it uniquely in accordance with Helmholtz theorem 
as the sum of two parts J = Jt + Jl where Jt is irrota- 
tional and Jl is solenoidal. ^"'■^^ Then it follows from the 
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continuity equation and the time-derivative of Q that 
Jl = —EodEc/dt. The result (a) may be reformulated 
in a way similar to l|44|l but then we only get the left im- 
plication (i.e."<J=") and the last term in (|43|l is then as- 
sumed to be linear in time. The leaking capacitor is one 
of a few examples below showing that quasistatics may 
apply exactly also with a nonlinear time-dependence of 
the current. 

Let us now prove H44() and start with the right im- 
plication. That Darwin is satisfied follows directly from 
the equivalence (|2U|I . Then both Ampere-Maxwell and 
Ampere-Darwin are satisfied implying dEp/dt — 0. 
From © it then follows that d^B/dt'^ = 0. From the 
second order time-derivative of Ampere-Darwin and H12|l 
we now get 

and (|43|) easily follows (including the conditions V • a = 

V - b = 0). Consider next the left implication in (|44|l . By 
substituting (|43|l in the Maxwell-Darwin equation we find 

V X B = /io(at + b) so that B is at most linear in t. From 
the time-derivative of © we then get dEp/dt = 0. In 
this case the Ampere-Darwin and the Ampere-Maxwell 
laws are the same and the solution to Darwin also solves 
Maxwell's equations. 

Above we found a necessary and sufficient condition 
(I44|l for the exact validity of Biot-Savart's law within 
Maxwell's equations. Thereby the result (a) of Griffiths 
and Heald was generalized. Let us now consider the cor- 
responding result for Coulomb's law and allow for a cur- 
rent of the form 

J(r,t)=b(r)-So^^|^ (46) 

where V • b = 0. We will prove that 

Maxwell + Coulomb EQS + Eq.^ (47) 

Let us start with the right implication in (|47ll . From 
Coulomb's law it follows the E-field is conservative so 
that EQS is obtained from Maxwell. From Faraday's 
law it then follows that dB/dt = and from the time- 
derivative of Ampere-Maxwell we now get (using E = 
Ec) 

and 14t)|) easily follows. To obtain the left implication 
of (|47|l we first substitute 146() in Ampere-Maxwell and 
obtain. V x B = ^oh. This equation together with 112|l 
results in a time-independent B-field and thus Faraday's 
law is satisfied. This is what is needed for the Maxwell 
equations to be valid for an EQS solution. The proof of 
()47|l is now completed and the result (b) is generalized. 

Let us finally consider the result (c) involving Ampere's 
law. This condition turns out to be not only sufficient 



but also necessary. Consider a current density of the form 

J(r,i) =a(r)t-f b(r) (49) 

where the vectorfields a and b satisfy V • a = V • b = 0. 
Then 

MaxweU 4- Ampere MQS + Eq.^ + Eq.Q 

(50) 

Consider first the right implication. Both Ampere's law 
and Ampere-Maxwell are valid so dE/dt — 0. By us- 
ing the time-derivative of Faraday's law we then find 
d^B/dt^ = so from the second order time-derivated 
Ampere's law we get d'^J/dt'^ = 0. Then (|^ follows if 
we also use that the current has zero divergence. The 
left implication in (|47|) is obtained by just reversing the 
procedure above. 

We will now consider a few examples involving the re- 
sults above. The first one will also be used in the next 
subsection. 

Example 1 Consider an (infinitely) long solenoid or a 
toroidal coil with time-varying current I{t). We assume 
the coils are winded so that we may neglect the axial cur- 
rent in the solenoid and the poloidal current in the toroid. 

Within quasistatic we get (exactly) zero B-field out- 
side the coils (from the Biot-Savart law and symmetry). 
The MQS and Darwin models coincide because there is 
no charge density and accordingly Ec = 0. There is 
a magnetically induced electric field E = E^ outside 
the coils in spite the fact that the B-field vanishes ex- 
actly. In general this is only an approximative solution 
to the Maxwell equations since the Ampere law (and of 
course the Ampere-Darwin law) rather than the Ampere- 
Maxwell law is satisfied outside the coils. An exact so- 
lution is obtained, as we have seen, when the current is 
linear in time. The electric field is then constant in time. 

The second example shows that quasistatics may agree 
exactly with Maxwell's equations even when he time- 
variation is quite arbirary. 

Example 2 Consider any given spherically symmetric 
charge and current density satisfying the continuity equa- 
tion 

By the spherical symmetry we get VxE = VxB = 0. 
The solution to EQS, Darwin and Maxwell's equations is 
E — Ec, B = and J — —eodEic/dt which is valid for 
any prescribed charge density p = p{r,t) in ©. The 
equivalences (|44|) and (|47|l apply with a = b = 0. A 
particular case is a leaking spherical capacitor. We as- 
sume that the medium in between the spherical plates 
has some conductance and does not violate the spheri- 
cal symmetry. In between the plates there is, during the 
slow discharge, zero charge density. The Ampere law is 
clearly not satisfied because B = while J ^ 0. 

The examples 3-5 below illustrates that in a linear, 
isotropic and homogeneous conductor the conduction 
current (assuming Ohm's law), in certain cases, does not 
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create any magnetic field. This corresponds to exact so- 
lutions of Maxwell's equations that also solves EQS and 
Darwin. 

Example 3 Let us now consider an example without 
symmetry in the charge distribution. Space is assumed to 
be linear, isotropic, conducting and homogeneous charac- 
terized by (eq, /j.0, tr). Assume that we at time t = know 
the charge density p(r, 0) = po (r)- 

From J = (tE in the continuity equation and by use of 
Gauss law we find 

p(r,i) =e-('^/^'°)Vo(r) (51) 

The corresponding solution to Maxwell's equation is 

E(r,t) = e"(''/"''^*Eo(r), B = (52) 

where Eg (r) is the Coulomb field due to the charge den- 
sity po (r). 

Example 4 The last example above may be modified so 
that we are dealing with a finite conductor surrounded 
by free space. Let the conductor contain all given charge 
Po (r) at time t = and let the surface of the conductor 
be an equipotential surface of Vq (r) defined by 

Inside the conductor both the charge distribution and 
the solution to Maxwell's equations is the same as in 
the previous example. Outside the fields are static with 
E = Eq (r) = —Wq and B = 0. The surface charge 
density on the conductor is 

as (r, t) = eo (l - e-^^/'°^') n (r) • Eq (r) (54) 

where n is the outward normal. 

Example 5 Related to the last example is the leaking 
capacitor of arbitrary geometry and with homogeneous, 
isotropic and weakly conducting material in between the 
two perfect conductors.^ 

In this case we must define Vq (r) by Laplace equation 
with appropiate boundary conditions on the two perfect 
conductors. The solution may be then expressed in terms 
of E = Eq (r) — — VVq as in the two last examples above. 

V. A QUASISTATIC PERSPECTIVE ON 
QUESTION #6 

A time-varying current of a long solenoid causes an 
induced electric field at the outside. But how can that 
be? There is no magnetic field at the outside! A.P French 
asked a similar question and it was addressed in several 
papersi^iSiLSiS The answers included two essential points 



• The vanishing of the B-field is a static phenomena 
while in the nonstatic case there is a small magnetic 
field outside the solenoid. 

• It is not good to consider the time-varying elec- 
tromagnetic field as a source for the induced elec- 
tric field, the source is rather the time-varying 
current 1^ 

However, even though both these statements are cor- 
rect in view of the Maxwell equations, there may remain 
some uncomfortable feelings. The way we in practice cal- 
culate the induced electric field from Faraday's law make 
it natural to think in terms of the time- varying B-field as 
a source and, even if there is a small B-field outside the 
solenoid, the main part of this source is well separated in 
space from the effect we consider (i.e. the induced out- 
side E-field). This apparent action at a distance seems 
to contradict the local action quality of Maxwell's equa- 
tions. But of course, a careful analysis shows that there 
is no true contradiction involved. 

Let us now consider this problem from a different point 
of view. If we use a quasistatic model in the calcula- 
tions it may be favorable to also think in terms of this 
model and not in terms of the full set of Maxwell's equa- 
tions. This is analogous to what we often do in other 
fields of physics. For example we use classical mechanics 
without worrying about quantum mechanics, Newtonian 
kinematics without referring to special relativity or New- 
tons law of gravitation without thinking about general 
relativity. It is in the solenoid example easy to calcu- 
late the magnetic and electric fields using quasistatics. 
One gets an exact solution within the quasistatic model 
and a good approximation to the true Maxwell solution. 
This quasistatic solution (of MQS and Darwin) has no 
magnetic field at all outside the solenoid. This solution 
is fundamentally inconsistent with Maxwell's equations 
since there is a time-varying E-field in free space where 
the magnetic field vanishes exactly. However, the so- 
lution is of course consistent with quasistatics. This is 
an indication of the quite different nature of quasistatic 
models and Maxwell's equations. The PDE's appearing 
in quasistatics are often elliptic or parabolic while the 
Maxwell equations are hyperbolic. The conceptual prob- 
lems have their roots in using quasistatic calculations but 
interpreting the result using Maxwell's equations. 

VI. ELECTROMAGNETIC SIMULATIONS 

A. The PC in basic courses 

Electromagnetic simulations are essential to electrical 
and electronic product designs in many industries. A 
broad range of important applications are within the 
quasistatic regime like motors, sensors, power genera- 
tors, transformer systems and Micro Electro Mechanical 
Systems (MEMS). Industrial and scientific applications 
may involve extensive calculations and accordingly the 
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need for much computational power. However, the usual 
PC has developed enormously and so has the software 
for simulations with increasingly user friendly interfaces. 
Today the PC is a sufficiently powerful tool for model- 
ing many electromagnetic phenomena within just seconds 
of computational time. This makes it potentially useful 
as an instrument for teaching basic electrodynamics. In 
the PC lab students may simulate various electromag- 
netic phenomena. Examples within statics or quasistat- 
ics include the edge effects in the parallel plate capacitor, 
distributed currents in conductors of various shapes, the 
Hall effect, magneto resistance, the appearance of non 
vanishing charge density in inhomogeneous conductors, 
the electric field outside and surface charge on a conduc- 
tor with currents, various objects placed in a given ex- 
ternal static or time- varying electromagnetic fields, eddy 
currents, inductive heating, magnetic diffusion, magnetic 
shielding and more. 



The use of electromagnetic simulations as an effec- 
tive pedagogical tool also presupposes some support from 
theory. One should in particular focus more on how 
Maxwell's equations are used to formulate well-posed 
PDE-problems for various physical situations. Most 
tcxtboks, at least the basic ones, are not much influenced 
by the appearance of computers. Their guiding principle 
is to find analytical solutions in simple and instructive 
cases. The underlying PDE-problem may not be needed 
explicitly in some examples like when integral laws are 
used in very symmetric cases or when tricks like mirror 
charges or mirror currents work. When the PDEs are ac- 
tually used, and this is much more the case in the more 
advanced texts, one finds solutions in terms of integrals 
and series by comparatively tedious calculations involv- 
ing, for example, special functions and variable separa- 
tion. 



With a numerical PDE-solvcr available the situation is 
somewhat different. The formulation of well posed PDE- 
problems for various phenomena is now motivated with- 
out the intention of finding analytical solutions. This 
is an easy part of the PDE theory for electromagnet- 
ics and it is instructive by providing increased insight 
into the structure of Maxwell's equations. It is in par- 
ticular important and sometimes straightforward to find 
how Maxwell's equations reduce thanks to various sym- 
metries. For example, the independence of one Cartesian 
coordinate or alternatively an axial symmetry results in 
a decoupling of Maxwell's equations. A simple inspec- 
tion of the equations written in component form reveals 
this structure. Such symmetries explain, for example, 
why the magnetic field of a toroidal coil only has an az- 
imuthal component or why there are TE and TM modes 
in planar wave-guides. In the next subsection we illustate 
this by formulating PDE problems for eddy currents. 



B. Equations for eddy currents 

Eddy currents and associated phenomena like induc- 
tive heating, magnetic shielding and magnetic diffusion 
are in practice very common since the basic ingredients 
are just a conductor and a time- varying electromagnetic 
field. Still eddy currents are not much discussed in the 
textbooks. An exception is the one by Smythe where a 
whole chapter is devoted to eddy currents and analytical 
solutions are given in terms of series and integralsj^i An 
interesting method of images may be used to calculate 
eddy currents in a thin conducting sheet. The theory 
was developed by Maxwell and reformulated in modern 
terms by Saslow.'^^ However, the analytical theory for 
eddy currents is still comparatively difficult. Electromag- 
netic modeling on a PC is today an attractive and easy 
way to include more about eddy currents in basic courses. 

Eddy current problems fall into two classes, steady- 
state and transient. In steady state analysis (also called 
time-harmonic analysis) we simply replace d /dt with the 
factor jbj and allow for complex valued fields. Maxwell 
equations becomes time-independent and much easier to 
solve. Physically we may obtain a steady-state condition 
sufficiently long time after the start of a time-harmonic 
source of the fields. The initial transient behavior is not 
considered in this analysis. The steady state Maxwell 
equations may be further simplified in quasistatic situa- 
tions, this may be useful for analytical theory but is not of 
much interest for electromagnetic simulations. However, 
in time-dependent (transient) analysis an initial value 
PDE must be solved. Then a quasistatic approximation 
may imply a major numerical simplification changing a 
hyperbolic PDE into an parabolic or elliptic one. 

Let us assume that the time-variation is slow enough 
so that quasistatic theory applies. But what quasistatic 
model should we use in the study of eddy currents; EQS, 
MQS or Darwin? Obviously not EQS since magnetic 
induction (Faraday's law) is outside the scope of that 
model. What about MQS? This seems to be the standard 
model for eddy currents calculation but it may in fact 
only be used in some particular casesii^i^L^ A quite com- 
mon misunderstanding concerning the quasistatic ap- 
proximation may partly explain the popularity of MQS.^j^ 

Let us now consider situations where the use of MQS 
may be justified. This happens in certain symmetric 
cases when MQS and Darwin are equivalent models. We 
consider below two examples of such symmetries. The 
first one include, as a special case, the situation when 
Maxwell's theory for eddy currents in thin conducting 
sheets appliesi^ 

Example 6 Consider a body in which the conductivity 
a = (7(z) only varies in the z-direction. In the xy— di- 
rections the body is homogeneous and of "infinite " extent. 
The electromagnetic field is created by the use of an ex- 
ternal current density 3'^ = Jxi^y 2/i •^)^ + ^' 2/' 
with vanishing divergence; V • J*^ = 0. 

In this example there will appear no charge density 
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and the Darwin model is equivalent to MQS. The reason 
is that electric fields will only appear in directions of 
constant conductivity (i.e. they have no z— component) 
and they cause no pile up of charge. The system may 
be described in terms of the potentials and Ay while 
V ^ A;, ^0. From equation ^ with J = J^+ai-W- 
dA I dt) we get the following equations for A^ and Ay 



dA, 



dt 



(55) 



These 3D equations determine the dynamics and may 
be used to model physical systems with the above sym- 
metry. 

A similar 2D axi-symmetric case may also be formu- 
lated. 

Example 7 We use cylindrical coordinates (r, 0, z) and 
a rotationally symmetric conducting body with conduc- 
tivity a — a{r,z). The external current density is 3"^ — 

.rit,r,z)l 

Also in this case there will appear no charge density. 
Some potentials vanishes, V = Ar ^ A^ ~ 0, and the 
system may be described in terms of the vector potential 
A = A{t,r,z)(l). This time we get the equation for the 
dynamics as 



dA Id, dA^ 



d^A 



dt 



-{r—) + 
r dr dr dz"^ 



(56) 



The equations may easily be studied by using a PC and 
some PDE-solver."^ For example, a straight copper wire 
above a copper plate is considered by Backstrom using 
equation (|55|l and the solver FlexPDE.'ii He also con- 
siders a similar example using H56(l for a circular copper 
wire above a circular copper plate. The transient fields 
appearing when a current is turned on in the wire is stud- 
ied (pages 75 and 99 in Backstrom's book). The same 
models are also easily studied using Comsols Mulphysics 
with the Electromagnetic Module. These problems are 
numerically simple and, of course, there are many other 
PDE solvers that may be used. Such numerical tools 
make it possible to consider many cases for which there 
are no analytical solutions. 

The examples above were carefully designed in order to 
avoid charge density. This is also true for all eddy current 
examples in Saslows' paper (see in particular Appendix 
B of that paper). 3^ Usually there will however appear 
time- varying surface charge on conductors with eddy cur- 
rents and, in case of an inhomogeneous conductor, there 
will also appear charge density inside the conductor. In 
order to describe the physics we then need to include 
the scalar potential in the analysis and MQS cannot be 
used An interesting possibility is then to use the Dar- 
win model. Such applications of Darwin's model have 
been suggested in connection with eddy currents in the 
human body caused by high voltage transmission lines. 4- 



The traditional use of Darwin's model is not eddy cur- 
rents but concerns charged particle beams and plasma 
simulations. 

Let us write equations for the Darwin model in term 
of potentials. The conducticity of a possibly inhomoge- 
neous conductor is cr = air) and that the time-dependent 
external current — J"^ {t, r) is prescribed The total 
current may be written J = -I- a{—W — dA/dt). We 
use, as always in the Darwin model, the Coulomb gauge 
V-A = 0. The dynamics of the potentials are determined 
by the continuity equation 



-eoTrVV - V • (aVV + adA/dt) 
ot 

and the Ampere-Darwin equation 



-V - J'' 



(57) 



d 



IIq{(tVV + adA/dt) + V^A : 



(58) 

These equations have a non-standard appearance by 
containing mixed time and spatial derivatives. They may 
however be solved using Comsol Multiphysics. Actu- 
ally there is a considerable freedom for the user to en- 
ter non-standard equations in this PDE-solver. This is 
achieved by making the equations available also on the 
"weak form" level which is the natural form for the fi- 
nite element method. The numerical solutions appears, 
at least qualitatively, to behave in an expected way for 
the few examples that we have considered. The practical 
usefulness of the Darwin model for eddy current calcula- 
tions remains however to be proven. 

A qualitatively new feature of Darwin's model (as op- 
posed to MQS or EQS) is the possibility of resonance. 
Since the Darwin model includes both capacitive and in- 
ductive phenomena it may in principle be used to model 
systems where the energy oscillates between the electric 
and magnetic fields. Consider, for example, a field theo- 
retic version of a LC-circuit. A suitable design for such 
an application may be a resonator in the form of a short 
coaxial cable where one end is short circuited by a metal 
plate and the other end is still electrically open but with 
increased capacitance created by the use of two close par- 
allel plates, connected to the inner and outer conduc- 
tor, respectively. This kind of resonator has applications 
in connection with electron beam devices at microwave 
frequencies (example 3.4.1 in the textbook of Haus and 
Melcher).^^ The resonance frequency of a LC-circuit is 
f — 1/ %/ LC but the frequency must be low enough not 
to violate the quasistatic assumption. In the above de- 
sign of a resonator the most essential method to achieve 
this (in the absence of dielectrics or magnetic materials) 
is to make the capacitance large and thus take the two 
capacitor plates very close to each other. 



VII. SUMMARY 

Maxwell's equations are fundamental for the descrip- 
tion of electromagnetic phenomena and valid for an enor- 
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mous range of spatial and temporal scales. The static 
limit of the theory is well defined and of course much 
easier. The electric and magnetic fields may in this limit 
be given by the laws of Coulomb and Biot-Savart. How- 
ever, as soon as there is any time-dependence wc should 
in principle use the full set of Maxwell's equations with 
all their complexities. Time-retardation is a fundamen- 
tally important but also a complicating feature. Using 
Maxwell's equations means in analytical theory that even 
if the effect is small it will not vanish and this makes the 
theory unnecessary complicated. In numerical analysis 
these effects, however small, may force us to use smaller 
timesteps (for numerical stability) and expensive calcu- 
lations. It is therefore useful to introduce quasistatic 
approximations in Maxwell's equations. The quasistatic 
models are also useful for a better understanding both 
of low frequency electrodynamics and for explaining the 
transition from statics to the general high frequency elec- 
trodynamics. This have been discussed in the present 
paper and below we list a few major points. 

(1) The quasistatic limit of Maxwell's equations is a 
kind of c — !■ CX3 limit obtained by neglecting time- 
retardation. The Darwin model is obtained if we 
use the Coulomb gauge. 

(2) The Darwin model involves both capacitive and in- 
ductive features but there is no radiation and the 
interactions are instantaneous. Poynting's theorem 
for this model shows that there is both electric and 
magnetic energy, but the electric energy only in- 



cludes the Coulomb part of electric field. 

(3) EQS and MQS may be considered as approxima- 
tions of the Darwin model. EQS includes capaci- 
tance but not inductance while MQS includes in- 
ductance but not capacitance. Poynting's theorems 
for these models show that there are only electric 
energy in EQS and only magnetic energy in MQS. 

(4) The law of Biot-Savart is valid within EQS, MQS 
and the Darwin model. 

(5) The law of Coulomb is of general significance for 
quasistatics (EQS, MQS and Darwin) as is obvious 
when we use the formulation in terms of force laws 
(subsection IIA) or use potentials (subsection IID). 

(6) The law of Ampere is not of general significance 
within quasistatics. It is valid only in MQS but 
not within EQS and Darwin^ 

(7) Galilei invariance of quasistatics is a somewhat del- 
icate issue. Galilei invariance structures may be 
defined for EQS and MQS but if we like the force 
to have the corresponding invariance then only the 
electric force is included in EQS and only the mag- 
netic force in MQSi2^ 

(8) Quasistatics has important applications in electro- 
magnetic modeling of transient phenomena. 
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